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We study the field equations of modified theories of gravity in which the lagrangian is a general 
function of the Ricci scalar and Ricci-squared terms in Palatini formalism. We show that the inde- 
pendent connection can be expressed as the Levi-Civita connection of an auxiliary metric which, in 
particular cases of interest, is related with the physical metric by means of a disformal transforma- 
tion. This relation between physical and auxiliary metric boils down to a conformal transformation 
in the case of f(R) theories. We also show with explicit models that the inclusion of Ricci squared 
terms in the action can impose upper bounds on the accessible values of pressure and density, which 
might have important consequences for the early time cosmology and black hole formation scenar- 
ios. Our results indicate that the phenomenology of f{R,R fl , v R ,lv ) theories is much richer than 
that of f(R) and f(R flv R llv ) theories and that they also share some similarities with Bekenstein's 
relativistic theory of MOND. 
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I. INTRODUCTION 



The existence of a cosmological constant is a chal- 
lenging problem for theoretical physics which has gained 
increasing attention since the discovery of the cosmic 
speedup [l| . If it is thought as vacuum energy, a naive 
quantum mechanical analysis raises the question of 
why it is so small. If quantum field theory in curved 
space-times is invoked, then the question turns into 
why it is so large 0. If it is seen as having a quantum 
gravitational origin, then the dimensionless combination 
A(GH/c 3 ) < 1CT 123 suggests that this is 123 orders of 
magnitude away from the right solution. Interpreting it 
as a new fundamental constant could imply the existence 
of new physics in much the same way as the Planck 
constant did. For these and many other reasons, a 
strictly constant cosmological term is undesirable. 

Different versions of dark energy sources have been 
proposed in the literature aimed at relaxing the con- 
dition of strict constancy of the so called vacuum 
energy. A scalar field with very low kinetic energy can 
closely mimic an effective cosmological constant, though 
away from the minimum it might exhibit characteristic 
features that could help distinguish it from a purely 
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constant term. In the arena of modified theories of 
gravity, there exists a family of theories in which a 
strictly constant cosmological term arises in certain 
circumstances. We refer to the so called f(R) theories in 
Palatini formalism, which are dynamically inequivalent 
to their (non- identical) twin brother f(R) in metric 
formalism. In these theories, the field equations for 
any lagrangian f(R) become exactly the same as those 
of General Relativity plus an effective cosmological 
constant whenever the trace of the energy momentum 
tensor is constant. When the trace is not constant, one 
gets modified dynamics induced by the time and/or 
spatial derivatives of T. It is this property of generating 
an effective cosmological constant that first attracted 
attention on these theories in relation with the cosmic 
speedup problem [3[. Any f(R) theory which at low 
cosmic densities (T — * 0) generates a cosmological 
constant of the right magnitude will produce the desired 
accelerating effect. However, the f(R) models with 
infrared corrections considered in the literature so far, 
besides producing late time cosmic speedup, also lead to 
catastrophic effects at microscopic scales. Such models 
induce severe instabilities in atoms due to the strong 
gravitational backreaction that occurs near the zeros of 
the atomic wave functions. In those regions the energy 
density drops virtually to zero and that excites the 
characteristic infrared scale of the gravitational theory, 
producing undesired effects which disintegrate the atom 
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The negative features of infrared corrected models, 
however, are not present in models with just high cur- 
vature corrections. Moreover, there is no fundamental 
reason that prevents an f(R) theory with high curvature 
corrections from producing an effective cosmological 
constant compatible with observations. We believe 
that that line of research has not yet been sufficiently 
explored. On the other hand, high curvature effects have 
already been observed in Palatini f(R) theories. In fact, 
it is possible to find an f(R) Palatini lagrangian which 
avoids the Big Bang singularity. In [5| it was shown that 
the effective dynamics of Loop Quantum Cosmology 
(LQC) can be exactly reproduced by an f(R) Palatini 
theory consisting on an infinite series in R/R c , where 
R c ~ Rpianck |5[- That result is important because it 
establishes links between an approach to quantum grav- 
ity based on non-perturbative Hamiltonian quantization 
techniques and a covariant action in Palatini formal- 
ism. It is well known that perturbative quantization 
techniques require the addition of quadratic curvature 
invariants R 2 , R IU) W iV , . . . (in metric formalism) to the 
Einstein-Hilbert lagrangian to keep divergences under 
control. However, it was not known what kind of 
effective actions, if any, could be associated with the 
non-perturbative loop quantization until the results of 
Q were published. 

Motivated by those results, one may wonder about the 
role that other curvature invariants besides R might play 
in the dynamics of Palatini theories. Do R^ V R^ V terms 
introduce new dynamical effects not present in f(R) 
theories? Another practical question is how such terms 
could affect or contribute to the form of the resulting 
effective cosmological constant in the vacuum limit of 
the theory. In other words, can we find new insights 
or mechanisms to obtain a cosmological constant of 
the right magnitude at low energies? To answer these 
questions, we are forced to go beyond the simple f(R) 
models. That is the goal of this paper. 

Palatini theories with Ricci squared terms have been 
already studied in the literature. In [6(, the cosmology 
of f(R) and f{R jJiV R IJ,v ) theories was considered in 
some detail. It was found that the scalar R^R^ V 
could be expressed as a function of the trace T of the 
matter energy-momentum tensor and that new dynamics 
arises in a way similar to that of f(R) theories, where 
R can also be expressed as a function of T. In Q, 
theories of the form R + f{R^ v R pv ) rather than simply 
/(Rfj.uR^) were considered, and the main focus was 
on the cosmology at the first-order perturbation level 
(density perturbation growth). In that work a 3 + 1 
decomposition was introduced in order to solve the 
connection equation, and the solution was given at first 
order in the metric variables. In this work we go a bit 
farther and consider Palatini theories in which the la- 
grangian is a generic function of the form f(R, R^R^). 
We study their field equations and focus on the basic 



manipulations that will allow us to write the theory in 
a form suitable for applications. To be precise, we show 
how to exactly solve the connection equation in terms of 
the metric and the matter sources (without introducing 
the 3+1 decomposition used in [7[), and discuss how the 
scalars R and R^R^ can be expressed as functions of 
the density and pressure. We illustrate this point with a 
family of exactly solvable models and point out that the 
value of R^vRv in such models is bounded from above, 
which also sets upper bounds on the energy density and 
pressure. We comment on the possible consequences 
that this aspect might have for the early time cosmology 
of such models. We also find that R and R^R^ 
are not just functions of the trace T, as it happens 
in the simpler f(R) and /(R^R^) theories. This 
fact implies that f(R,R^R^ v ) theories have a much 
richer phenomenology during the early universe and, in 
particular, during the radiation dominated era. This is 
so because in that era, which is characterized by T = 0, 
the Palatini corrections of f(R) and f{R l i U R tl ' J ) theories 
boil down to just a cosmological constant, whereas in 
the f(R, R^vR^) case the modified dynamics is more 
involved. 

The paper is organized as follows. In section [TTl we 
illustrate how to solve for the connection in general 
f(R,R pv R pv ) theories. We start with the f(R) case to 
describe the basic steps of the algorithm and then focus 
on the general f(R, R^R^) case applied to a perfect 
fluid. Then we write the field equations of the metric in 
a form suitable for applications and discuss the vacuum 
limit. In section [TVl we provide a family of models for 
which R and R fliy R flL ' can be explicitly solved in terms of 
the density p and pressure EE, and show that those quan- 
tities must be bounded from above for consistency of the 
equations. We conclude with a summary and discussions. 

II. SOLVING FOR THE CONNECTION 

The action that defines a Palatini f(R, R^ V R^ V ) theory 
is as follows 

S[g,T,ip m ] = j ' d i x y ^gf(R,Q) + S m [g,ip m ] (1) 

where g a p represents the space-time metric, FjjL is the 
connection (which is independent of the metric), ip m rep- 
resents the matter fields, R = g^R^, and Q — R^R^ . 
The Ricci tensor R^ v — R tip ^ p is defined in terms of the 
connection as follows 

iV(r) = -dA + dxf% + r^r p pX - r* p r£ A (2) 

This definition follows from the relation [V a ,V p ]wp = 
Rapf3 S ws and is valid for any derivative operator V Q 
In order to obtain the field equations, we must vary the 
action with respect to the various fields present in it. 
For completeness, in Appendix we detail the process 
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of variation to get the field equations. Variation with 
respect to metric and connection lead to the following 
equations 



(3) 
(4) 



were we have used the short-hand notation fn = <9r/, 
and /q = <9q/. We now focus on working out a solution 
for g]). 

At first sight, equation ([¥]) seems a highly non-trivial 
equation because of the terms fn and /q-R^". In fact, 
since R and are functions of the connection and its 
first derivatives, (j4|) can be seen as a non-linear, second- 
order equation for the unknown connection. However, 
there exist algebraic relations between i?, R pv and the 
energy-momentum tensor of the matter that will make 
the problem easier. To illustrate the general algorithm 
that solves (UJ), we will first focus on the f(R) case, which 
is just the case /q = 0. 



A. Solving for the connection in /(-R). 

In this case, the connection equation becomes 

V/j [V^fng^] = . (5) 
From the metric field equation ([3]), we find 

IrR^ ~ = k 2 T m „ . (6) 

The contraction of this equation with the metric yields 

f R R-2f = n 2 T. (7) 

This expression is an algebraic equation that general- 
izes the relation R = —k 2 T of GR to arbitrary la- 
grangian f(R). This means that given a lagrangian, say 
f(R) = R + ai? 3 , the above equation can be solved as 
R = 1Z(T). In this particular example, the equation is 
just -R + aR 3 = k 2 T. Note that 1Z(T) is now a function 
of the matter and, for this reason, f(lZ) and f R (TV) are 
also functions of the matter. According to this, the con- 
nection equation ([5]) can be seen as a first order equation 
for the connection that depends on the metric g a fj and 
the matter: 



^-gf R {n{T))g^} 







(8) 



We now discuss how to solve this equation for the connec- 
tion. In GR, this equation is simply V/3 [>/—gg ,iU ] = 0, 
and the solution @ is given by the Christoffel symbols 
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up 



(9) 



In order to find a solution for ([5]), we can do the follow- 
ing. We assume that there exists a metric h pv such that 



the connection that solves ((8|) is the Levi-Civita connec- 
tion of h pu . This means that the metric h pl , satisfies 
V/3 [%/— hh* 1 "] = 0. In other words, our ansatz satisfies 



-gfW 



(10) 



It is useful to rewrite this equation using matrix notation 
to suppress indices: h pv — > h, h^ u — » h~ x . We then have 

V^hh- 1 = J—gfng- 1 (11) 

If we compute the determinant of the left and the right 
hand sides, we find 



gfk (12) 



Introducing this result back into (fTTj) . we find 

h^ 1 = g^ 1 1 'fn, or equivalently, = g^ / fn and 

h pv = fng^v We thus see that the metric h pv that 
defines our independent connection is conformally 
related to the spacetime metric g pv . 

At this point, we must make some clarifications about 
our notation. First of all, since we are dealing with two 
metrics, we can also construct two different Ricci tensors 
and two different scalar curvatures. The curvature 
of the metric g pv is defined as R(g) = g^Rpv{g) 1 
where R pv (g) is the Ricci tensor constructed from the 
Levi-Civita connection of g pv . Analogously, we can 
define the scalar curvature of the metric h pv as the 
contraction R(h) = h^ v R^(h), which is equivalent to 
bf"* i? M „(r), where T is the Levi-Civita connection of 
h pu . Now, in the above equations we have been dealing 
with something different from R(g) and R(h). What we 
called R or 1Z(T) is a hybrid object, it is the contraction 
g^ RpuiT) , which is defined in terms of the independent 
connection via R^iT) = Rp V (h) and the spacetime 
metric g^ v . In all our previous formulas, we dealt with 
the object 1Z, which is a function of the matter via the 
trace T. Whenever we write TZ it should be understood 
that it is a function of the matter. In any other case, we 
will write either R(g) or R(h). 

For completeness, let us now focus on the equation for 
the metric. In order to write a second-order equation for 
the metric, we must use the solution that we just found 
for the connection and introduce the result in ©, which 
is written in terms of R^ u (h) and g pv . Since h and g are 
conformally related, it is straightforward to relate R^(h) 
with R^ v (g) (see Appendix D of @). The result is that 
([6]) becomes 

G>„(flO = -r-Tnv - — ^ — —Qiiv + T - (VjuViz/r — p^D/r) 

JR ^JR JR 



3 / 1 

[dpfRdvfR - ~g^{df R f 



(13) 



where now V M is the usual covariant derivative of g^ v . 
The lesson to learn here is that the right hand side of 
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(I13p behaves like a modified energy-momentum tensor in 
which the trace T plays a role non-existing in GR via 
the terms 1Z, f(TZ) and the various derivatives of Jr(JV). 
When T is constant, all d^f terms vanish and the equa- 
tions boild down to those of GR plus an effective cosmo- 
logical constant given by evaluated at constant 

T. 



B. Solving for the connection in general f(R,Q) 

We now face the problem of solving the constraint 
equation In this case, the metric variation leads to 

f R R^ - ~g„v + 2f Q R^K = k 2 t^ (14) 

In order to solve ((4]) we will follow the same steps as 
above but adding the complexity inherent to the new 
problem. 



• Step 1. We need to find an algebraic relation be- 
tween R and R^ v with the matter sources. This will 
allow us to re-interpret as a first order equation 
for the connection. To proceed, we first define the 
matrix P, whose components arc P^ = Rfj, a g av , 
which allows us to express (|14p as 

/«P/ - £ + VqP^Pc? = K a V . (15) 
In matrix notation, this equation reads 

2/qP 2 + SrP - tl = k 2 T , (16) 

where T is the matrix representation of . Note 
that R and Q are the trace of P and the trace of P 2 
respectively. Solving this equation will thus lead to 
a relation of the form P = P(T), which is analogous 
to the solution TZ(T) of the f(R) case. For the 
moment we will not care about the particular form 
of this solution, which depends on the particular 
model chosen, and will just assume that it exists. 
This is what we need to reinterpret ^ as a first 
order equation for the connection. The solution 
should thus depend on the metric and the T^ v 
of the matter. 

• Step 2. We now propose an ansatz of the same 
form as above. We look for a metric h such that 
V/3 W~ hh* 1 "] = 0. This guarantees that the con- 
nection will be the Levi-Civita connection of h. Us- 
ing matrix notation, we have 



-hhr 



-gg 



- l (f R I + 2f Q P 



(17) 



We now compute the determinant of the left 
and the right hand sides, which give h = 



g dot + 2/qP \ . Once we know the explicit 

expression for P we will be able to compute this 
determinant. In any case, we have the formal ex- 
pression 



(18) 



Vdett 

where we have defined E = (JrJ + . Taking 

the inverse of the above matrix, we find 

h = (VdctE) ±- x g . (19) 

We have thus shown that the connection of f(R, Q) the- 
ories can be explicitly solved in terms of the physical 
metric g M „ and the matter sources. To proceed further, 
we need to consider particular choices for T^ v and/or 
particular models. This is the task of the next sections. 

C. f{R,Q) with a perfect fluid 

We now consider the case of matter described by = 
n<Sj t 1 ' + (/9+n)w /J w l/ , where II represents pressure to avoid 
missunderstanding with the matrix P. The first thing 
that we need to do is to compute the matrix S and its 
inverse, because that will allow us to compute h and h^ 1 . 
To find E, we first need to find P, which is a solution of 
([To]) . That equation can be rewritten as follows: 



fn A ( ,.2-n , f , /fl 1 / . n ,,/< 



4/q y V 2 

which can formally be denoted as 

2/qM 2 = al + f3u^u v , 
so that we can look for a solution of the form 



(20) 
(21) 

(22) 



where fg > has been assumed. We can now take the 
square of this matrix to find the relation between (A, a) 
and (a, ft): 



(-X±y/X 2 +pn u ) 



(23) 



where n u = u^u^ is the norm of the vector u^, which for 
the perfect fluid is just n u = — 1. If /q < 0, then the same 
manipulations hold up to the redefinitions /q — > — |/q|, 
a — ► —a, and j3 — > —(3. The matrix E can then be 
written as 



y v 



(24) 
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where we have defined 



A 2 = ,[2j7>o 



(25) 
(26) 



The determinant of a matrix of the form (|24[) can be 
computed straightforwardly using the definition det M = 
e abcd M a M 1 b M 2 c M 3 d and leads to 



detS = A?(Ai + n u A 2 ) 
The inverse (E -1 )^ has the form 

A 2 1 



-1^ 
Ai " 



Ai (Ai + A 2 n u 



-u^u 



(27) 



(28) 



We are thus ready to express the metric in terms of known 
quantities. 



hfxv ^ 



1 
9. 



<r + 



(29) 
(30) 



where = x/detE/Ai. It is worth noting that in the 
limit /g -> we find 91 —> fa and cr — ► when the plus 
sign in front of the square root of (|2"3")l is taken, which 
recovers the equations of f{R) theories. A relation be- 
tween two metrics of the form found in (|29p and (|3T))) . 
namely h^ u = Clig^u — Cl2U^u v , involving a vector field 
and two independent functions f2i and 9l 2 , is known in 
the literature as "disformal transformation". A relation 
of this form between two metrics has already been studied 
in inflationary models characterized by non-minimal cou- 
plings between curvature and a scalar field [To| , and nat- 
urally arises in Bekenstein's relativistic theory of MOND 
PH and in previous versions of it. In the MOND the- 
ory, the vector u p is an independent dynamical vector 
field and the functions in front of it and in front of g^ v 
depend on another dynamical scalar field. In the theory 
described here, on the contrary, the metric tensor is the 
only dynamical field of the gravitational sector. 



In the case of a perfect fluid, n v 
side of this equation becomes 



T, the right hand 



(4 + « 2 n) A lK 2 ( P + n)-A 2 (| 



2 n 



Ai ' 9 ^' Ai(Ai - A 2 ) 

The field equations for the metric g^ v can now be 
obtained from (|32|) by expressing the connection Y(h) in 
terms of g^ v and the matter using (f2"9")l and (|3"0"|) . In the 
case of f(R) theories, this process leads to (fl~3|) . (Recall 
that the f{R) field equations are recovered in the limit 
fg^O when the plus sign in front of the square root of 
(f2"3"| is taken.) In the case of f(R, Q) theories, we find it 
more convenient to work directly wit h (I32[) . Examples 
of applications will be given elsewhere [12j |. 

Let us now consider the effective cosmological constant 
A e ff that follows from a generic f(R,Q) theory in vac- 
uum. With the same choice of sign that recovers the right 
f(R) limit, (|32|) in vacuum becomes 



Rfi (r)-uac 



J_ 

2A X 



(34) 



p,n=o 



Since in vacuum the relation between and g^ v is 
just a constant conformal factor, h^ v = 9l\ pU=0 g^ v 
(with 9l\ pJI=0 = Ai| pn=0 ), it follows that R^ v {T)vac = 
R^v{g)vac, which implies that 



Rfj.u{g)vac — 2^ 



p,n=o 



9^ 



= A 



eff9nv 



(35) 



This equation indicates that A e f f is the result of evalu- 
ating the lagrangian / in vacuum, up to a factor l/2Ai. 
Now, from the trace of (|16p we see that in vacuum 
f/2 = {Rfn + 2Q/q)/4. Since we assume lagrangians / 
which recover GR at low curvatures, which implies that 
fn — > 1 in vacuum, it follows that to get a non-zero 
A. e ff we must have R vac 7^ 0, or QJq 7^ or a combi- 
nation of both. In any case, we see that the inclusion of 
Q-dependent terms in the lagrangian provides new mech- 
anisms besides those present in f(R) theories to generate 
a non-zero cosmological constant. 



III. FIELD EQUATIONS AND A eff 

To find the field equations for the metric in a way that 
reminds those of GR, we can take (| L 5 [) and put it as 

P^ a v = k 2 T^ + f -5^ . (31) 

Using the inverse of E and the definition of P M Q = 
R^ v g va , we can express these equations as 

R^(T) = (K 2 T^ a + ^6^)^-%%^ = r MV . (32) 



IV. SOLVING FOR R(p,U) AND Q(p,n). 

In order to study physical predictions of the theories 
considered in this work, we must specify particular mod- 
els. Once a model is chosen, the problem of finding 
R = lZ(p,H) and Q = Q{p, II) is difficult in general be- 
cause it can involve non-linear equations with multiple 
solutions. In principle, one only needs two independent 
equations relating R, Q with p, II but it will not be obvi- 
ous always how to construct two such equations to make 
the problem analytically tractable. Here we will impose 
some restrictions on the lagrangian to simplify and illus- 
trate the problem. We will consider f(R, Q) functions of 
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the form f(R, Q) = f(R) 



R U „R" 



where Rp is a con- 



stant of the order of the Planck curvature. The reason 
for this choice becomes apparent when we take the trace 
of equation (fTHJ) 



2Q/q + Rf R - 2/ = n 2 T 



(36) 



This equation relates R and Q with T, but in the partic- 
ular case of f(R, Q) — f(R) + R^R^ V j Rp, it becomes 



Rf R - 2/ = k 2 T 



(37) 



which is exactly the same expression as in /(i?) theories 
and implies that R = TZ(T). There still remains to obtain 
Q = Q(p, II). This function can be found taking the trace 
of I 



( R + ^ ) = 3A + 1 



which can be cast as 



3A 



= A 2 



(38) 



(39) 



After a bit of algebra we find that 



3 [R + 




(40) 

For our choice of f(R, Q) we see that the /q term ap- 
pearing on the right hand side of (|40|) is just the con- 
stant Jq = 1/Rp. We can then use the definition of A 2 
[see (pO]l. (j2i"j) , and ([22]>] to solve for Q as a function of 
lZ(T),p, and n. The solution is 



Q 

2R P 




R 
Rp~ 



± 



n) 



Rp 



(41) 



The correct limit when Rp — > oo is obtained for the mi- 
nus sign in front of the square root. Had we chosen Rp 
negative, i.e., Rp — > —Rp, (|4"Tj) would still hold but the 
positive sign in front of the square root should be needed 
to obtain the right limit. This can be easily verified for 
particular choices of /. 

It is remarkable that (|4"Tj) with Rp > implies that Q 
is bounded in the real line. Consistency requires that the 
term under the square root be non negative, i.e., 



R-p +fR ) — ; rp~ 



n) 



> o 



(42) 



which places bounds on the accessible values of p and 
II. To illustrate this point, let us consider the family of 
lagrangians f(R, Q) — R+aR 2 /Rp + Q/Rp, for which R 



turns out to go exactly like in GR, 1Z — —k 2 T. If in this 
example we choose a = —1/2, i.e., f(R) = R— R 2 / (2Rp), 
then ([4T|) becomes 



3R 2 P 



1 - 



2K 2 (p + n) 2« 4 ( / 9-3n) 2 



Rp 



3R 2 P 



I 4k 2 (p + TI) 
Rp 



(43) 



At low energies, this expression recovers the GR limit 
3(P + pf 15{P + p) 4 



(3P 2 



P 



2R P 



4R P < 



(44) 



However, positivity of the argument in the square root 
of (|43]) implies that at high energies 



^(P + U)<^f, 



(45) 



which clearly shows that the combination p + II is 
bounded from above. Taking an equation of state of the 
form II = wp, we find that the maximum value of Q 
occurs at n 2 p max = Rp/{^ + iu>) and is given by 



3i?p 

"TfT 



i 



(1 - 3w) 2 
12(1 + w) 2 



(46) 



It is important to note that the upper bound on p and 
II is a consequence of the field equations, which have 
been used to obtain (|4"Tj) . and is independent of the sym- 
metries or equation of state that our particular problem 
might have. As a consequence, we may expect to ob- 
serve important new effects in critical scenarios such as 
in stellar collapse processes or in the very early Universe, 
where singularities are unavoidable in the context of GR. 
We would like to remark that, within the metric formal- 
ism, theories with bounded curvature scalars were pro- 
posed time ago in [1 31 ] to construct non-singular isotropic 
universes. In such constructions, the Big Bang singular- 
ity is cured by the self-interactions introduced by two 
scalar fields, whose role was to place bounds on the up- 
per limits accessible to certain curvature scalars. In the 
theory discussed here, no new degrees of freedom have 
been introduced. It is the Palatini dynamics that con- 
straints the physical range of p and II and, therefore, 
suggests that all curvature invariants will be bounded. 
This contrasts with the generic behavior of the models 
R+ (aR^R^ + bR 2 )/Rp in metric formalism, where all 
solutions which at late times recover a FLRW Universe 
start with a Big Bang singularity and, therefore, with 
divergent scalars R and Q |14l |. 



V. SUMMARY AND CONCLUSIONS 

In this work we have studied the field equations of 
Palatini theories of gravity in which the lagrangian is a 
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function of the form f(R, i? M „i? MI/ ). We have shown that 
the independent connection can be expressed as the Levi- 
Civita connection of an auxiliary metric which is related 
to the physical metric g^ and the energy-momentum 
tensor by means of a non-standard transformation, which 
becomes disformal when matter is described as a per- 
fect fluid and boils down to conformal when the R^R^ U 
dependence disappears. The emergence of two metrics 
related by a disformal transformation, a basic require- 
ment of relativistic MOND theories to properly account 
for gravitational lensing, could make these theories inter- 
esting for the consideration of dark-matter-related prob- 
lems. 

We have also shown that in Palatini f(R, R tlv R ,iV ) the- 
ories the scalars R and R^ V R^ U can in general be writ- 
ten as functions of p and II but not necessarily via the 
trace T. As a result, the phenomenology of these theo- 
ries is much richer than that of the individual f(R) or 
f(R^ u R^ u ) theories. Furthermore, for some simple mod- 
els, we have shown explicitly that the scalar R^R^" sets 
bounds on the physically accessible range of p and II, 
which suggests that scenarios such as the very early Uni- 
verse and the last stages of stellar collapse could be se- 
riously affected by the new dynamics possibly leading to 
singularity resolution. 

The results obtained in this work open new avenues of 
research in the context of the very early Universe, the ra- 
diation dominated epoch, and the accelerating Universe, 
with new mechanisms to generate an effective cosmologi- 
cal constant, within a framework that does not introduce 
new degrees of freedom and, therefore, is closer to GR 
than other type of modified theories of gravity or dark 
energy models. 
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Appendix A: VARIATION OF THE ACTION 



and the subindex in f R and Jq denotes partial derivative. 
Since R — g^R^ and Q — g^ a g vl3 R^R a (3, it is easy to 
see that 

SR = 5{g^R pv )=R llv 8gi iV + g llv 5R llv (A3) 

= 2R tia R^ v + 2R^SR lil/ (A4) 
Inserting these results in SS we find 

SS = i / ^Xyf^[{f R R pv - f - 9llv + 2f Q R iia RZ)5g^ 



(f R g^ + 2f Q RnSR^} 



(A5) 



The next step requires to express SR^ in terms of <£T" . 
This can be done using the so called Palatini identity 



(A6) 



We now manipulate the SR^ term. That contribution is 
of the form 



M = J (Pxyf^gWdR^ 



(A7) 



where, in our case, A^ = (f R g» v + 2f Q R» v ). Using the 
Palatini identity, we get 



M 



dV^A"" [-V„ (STl) + V A (51%,)] (A8) 



Using integration by parts and rearranging indices, we 
find 

M = J d 4 x{v x y=g [A^ST X V - A Xv 5T^)] + 

V4^(A^^-^)]5rt}(A9) 

The first term in brackets is a total derivative and can be 
discarded. The second term is the one we need. Putting 
this back into 5S we end up with 

SS = ^ J d 4 x^{(f R R^ - f - 9llv +2f Q R lia R«)8g» v 

+V M [/=g(A^ - 5»A X »)] «l} (A10) 

Knowing that the matter action gives SS m = 
— h J d^XyJ—gT^Sg^, the field equations can be writ- 
ten as follows 



The variation of the action Q]) can be expressed as 



-tg^Sg^ + 5f(R,Q) 



where Sf(R, Q) represents 

5f(R,Q) = f R SR + f Q 6Q 



, (Al) 



(A2) 



f 

f R Rnv — -^g^v + 2fQR IMX R™ — k T M „ (All) 



~g (k^5 X p ~ ^A Aiy ) 



= 



(A12) 



Note that the second equation is equated to zero because 
matter is not coupled to the independent connection (<5r 
does not appear in 5S m ). Equation (|A"l2l) can be further 
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simplified if one notices that when A = j3 then the equa- 
tion is identically zero. Taking X =/= (3 then it boils down 
t° [ % / :r gA Al/ ] = 0, which is explicitly given by 

[V^ (finT + 2f Q Rn] = (A13) 



It is straightforward to verify that when / = f(R), no 
explicit dependence on Q, all the above equations reduce 
to the case of Palatini f(R). 
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